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ABSTRACT 


Heat-transfer measurements were made for ae single, 
water-heated tube in a pool of R-114 to simulate operating 
conditions of awater chiller. Data were obtained for a 
smooth copper tube, and for two commercially available 
tubes: a spirally roped copper-nickel tube with a porous 
coating; anda copper tube with a structured outer surface 
and a multiple-start helical ridged inner surface. 
Measurements were made for refrigerant-oil mixtures at oil 
concentrations from 0 to 6 mass’ percent with a boiling pool 
temperature of 13.8 nGe Results for the two enhanced tubes 
with and without oil are compared to the smooth tube data. 
Enhancement factors for the overall heat-transfer coeffi- 
cient were 4.0 and 3.6 for the structured surface and 
porous-coated tubes, respectively, in pure refrigerant and 
at a water velocity of 2 m/s. For these same conditions the 
enhancement factors for the outside heat-transfer coeffi- 
cient were 14.6 and 6.4 for the porous-coated and structured 


surface tubes, respectively. 
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I. INTRODUCTION 


Recent developments in boiling surfaces have shown 
considerable enhancements in heat transfer performance. The 
worldwide literature on enhanced heat transfer contains over 
3000 published technical papers and reports [Ref. 1: p. 81]. 
This increased interest in heat transfer augmentation is the 
result of incentives for energy and material savings. One 
effective way to improve the heat transfer is by uSing 
passive augmentation. 

Passive augmentation uses fine-scale alteration of 
surfaces, both external and internal [Ref. 1: p. 82]. These 
surfaces may consist of either an applied porous coating or 
fins deformed in various ways to provide a large number of 
reentrant cavities on external surfaces. Methods of 
internal enhancement include insert devices, forged fins, or 
deformation of the surface (i.e., Spirally roped or corru- 
gated tube). The use of such surfaces can lead to consider- 
able reduction in the size and weight of heat-transfer 
equipment. The reductions may lead to smaller capital costs 
or operating costs or both. 

Of particular interest, in naval applications, is the 
reduction in the size of water chillers in refrigeration 
systems. While a number of investigations are currently in 
progress, the theoretical treatment of the boiling perform- 
ance of various tubes is almost impossible owing to the very 
complicated mechanisms involved. The degree of difficulty 


increases with the presence of oil in the refrigerant 


ieui1d. In general, refrigeration systems with oil- 
lubricated, hermetically sealed compressors contain small 
mass percents of oil in the refrigerant liquid. Therefore, 


reliable data covering a wide variety of operating condi- 


moms for various refrigerants and different boiling 


ie 


surfaces is essential to further develop more compact and 
efficient evaporators. 

Most of the experiments on enhanced boiling surfaces 
reported inthe literature have used electrically heated 
tubes. The difficulties associated with these tubes, espe- 
clally during the instrumentation stage, raise questions as 
to the reliability of the data as described by Wanniarchchi 
et al [Ref. 2: p. 14]. Despite the precautions taken to 
minimize contact resistance, it may still be present between 
the thermocouple locations and the outer boiling surface. 
Another area of doubt, described by Wanniarchchi et al. 
[Ref. 2: p. 14], is the nonuniform heat fluxes generated by 
commercially available electrical heaters. Additionally, in 
order to use an electrical heated tube, the inside surface 
of the tube must be smooth, or any internal enhancement must 
be bored out. This is disadvantageous as designers would be 
interested in knowing the overall heat-transfer coefficient. 
Most industrial experiments use 2-to 3-meter-long tubes 
with warm water flowing inside them. While the information 
generated from these experiments closely approximate actual 
evaporators, compared to electrical heated tubes, the large 
water temperature drop from inlet to exit (up to 5 K) may 
make it difficult to study the details of the beau 
process. 


The type of refrigerant selected will also influence the 


design of evaporators. R-114, a moderate-pressure refrige- 
rant, is receiving more attention, in particular for naval 
applications. Advantages to using R-114 include: (a) it 


belongs to the refrigerant group with the lowest toxicity, 
(b) it is very stable with temperature, and (c) it haswa 
fairly large value for the emerey UtCrancten per unit volume 
Of vepor (leah, 2 eine /Vg) (Ref. 33) pp ee ee The 
advantages of the first two items are clear. The third item 


is advantageous since a higher value for E, means’ lower 


2 


pressure drop through the refrigerant piping for a specified 
heat duty. As the normal boiling point of the type of 
refrigerant increases, the E, value increases. R-ll, R-113, 
R-114, and R-22 have E, values in increasing order. R-114 
may be preferable to R-ll and R-113 due to the higher E, 
value and preferable to R-22 due to the ability to use 
lighter components in the refrigeration loop. 

Based on the above discussion, the major objective of 
the present investigation was to study the boiling perform- 
ance of two commercially available tubes, a porous-coated 
tube and a mechanically structured surface, in comparison | 
with a smooth tube for the following conditions: (a) tubes 
are water-heated, (b) boiling fluid is R-114 with 0, 1, 2, 
and 6 percent by mass oil, and (c) boiling temperature is 
Jams °C. 


Is 


Il. BACKGROUND 


Ave NUGEEATE SB OREING 

When the temperature of a boiling surface exceeds the 
Saturation temperature of a fluid by a few degrees, nucleate 
boiling occurs. The difference between the boiling surface 
and saturation temperature is the amount of wall superheat. 
Dougherty [Ref. 4: p. 175] defined pool-boiling as vaporiza- 
tion occurring under the following conditions: (1) Game 
depth >> the bubble diameter, (2) there is negligible effect 
on heat transfer due to low or no externally imposed 
velocity of the fluid, are (3) the bubbles move away from 
the boiling surface due to a force field. With the above 


conditions, the heat flow can be expressed as: 


ae wee 
h mt Se ( 22 IL 


A description of the boiling mechanism is given by 
Chongrungreong and Sauer [Ref. 5: p. 701]. A thin layer of 
Superheated liquid is formed adjacent to the boiling 
surface. Bubbles nucleate in this thin layer and grow from 
preferred spots on the boiling surface. It 1S assumed tae 
the primary resistance to heat transfer iS Within theweame 
liquid layer. The height of the liquid in the “boris 
container is not a primary variable, but the liquid prose. 
ties should be controlled. However, large flooded evapora- 
tors may experience a significant "submergence effect" from 


the liquid head. 


3 EXTERNAL SURFACE ENANCEMB NT 
For smooth heating surfaces, bubbles nucleate at various 
Scratches and cavities on the surface. Fujiilt stated times 


number of active Sites inereaseS aS the heat flux ine. a ee 
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[Ref. 6: p. 48]. Webb [Ref. 7: p- 46] stated that the 
ability to enhance the nucleate boiling coefficient by 
applying some type of roughness has been known’ for over 50 
years. 

Presently, there are two major types of commercially 
developed roughnesses: CMyieporomsscoating and (2) various 
fins and surfaces with reentrant cavities [Ref. 8: p. 24]. 
The first type of roughness is a Sintered particle coating, 
usually copper or aluminum. Figure 2.1 shows a schematic and 
a photomicrograph of a porous coating. The small reentrant 
cavities are interconnected by substrate tunnels. The 
surface reduces the required superheat for vapor generation 
by entrapping a high density of relatively large vapor 
nuclei in the cavities contained within the porous coating 
[Ref. 8: p. 24]. As described by Webb [Ref. 7: pp. 58], 
researchers found that a critical pore size not a particle 
Size governs the number of nucleation sites. Large pores are 
required for liquids with high surface tension and high 
thermal conductivity, while small pores work best for 
liquids with low surface tension and low thermal  conduc- 
tivity (i.e., refrigerants). Czikk and O'Neill [Ref. 9: pp. 
53] developed correlations to include the effects of 
different pore sizes of porousS-coatings. They concluded that 
there were two resistances to heat transfer for these coat- 
ings: a nucleation superheat related to bubble diameter; and 
a conduction superheat controlled by the liquid film sepa- 
rating the vapor bubbles from the metal matrix. Carnavos 
[Ref. 10: pp. 106-108} found that a porous coating resulted 
in a /00-800 percent better heat-transfer coefficient than a 
Smoetn Surface in a pool of R-ll. Czikk [Ref. 11: p. 98] 
also found increased performance of the porous-coated tubes 
used in his 20-ton water-chiller experiment. 

Structured fins or surfaces form reentrant cavities of 


various geometries. Figure 2.2 is shows an example of a 
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manufactured reentrant cavity. These reentrant cavities act 
as stable nucleation sites, which enhances the heat 
transfer. The nucleation site, to remain active, is depen- 


dent on the mouth diameter falling within a critical range 


and shape with a maximum reentrant angle. This range is a 
function of the fiuid = propereies- The cost of these 
surfaces, presently, are not Significantly higher than the 


cost for smooth tubes and they dramatically improve boiling 
performance compared to the smooth tubes. 

The mechanisms by which the reentrant cavities operate 
are complicated. Previous reSearch concluded that the cembas 
nation of bubble evaporation, thermal boundary-layer strip- 
ping and bubble agitation controlled the heat transfer from 
smooth surbtaces li Ret. 12-5 pa) ale zae In contrast, difftereme 
experiments conducted by Arshad [Ref. 12: p. 192] and Arai 
(Ret. 133" Yor 37} show that thin film evaporation in the 
reentrant cavity controlled the heat-transfer mechanism. 
Bubbles were formed by vapor exiting cavities as the liquid- - 
vapor interface of the thin film contacted the cavity 
surface.. Surface tension holds most of the liquid on the 
cavity walls. Ayub [Ref. 14: p. 64] showed this similar 
"thermosiphon mechanism" with his experiments using enhanced 
Surfaces. 

Nakayama [Ref. 15: _ p. 37] gives a fairly detailed 
description of a physical model undergoing this’ bubble 
growth mechanism. Figure 2.3 shows the three major phases of 
bubble growth. Phase I consists of a pressure buildup in 
the tunnel by evaporation of liquid held in corners’) and 
continues until the meniscus reaches a hemispherical shape 
of radius r, = dg/2 (where d, is the mouth diameter of (ie 
reentrant cavity). In Phase II, the meniscus grows faster at 
Some pores than at others. The Phase I pressure buildup is 
reduced as the vapor enters the growing bubbles. The 


meniscus at inactive pores does not grow due to this vapor 
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pressure reduction. Initially, the bubble expands under the 
high internal pressure, while lacer expansion is controlled 
by the receding liquid inertia around the bubble. PHaseuia 
is termed the liquid intake phase with liquid flowing aaa 
the tunnel through inactive pores. This flow is over a short 
interval of pressure depression which occurs as the pressure 
of the bubble and tunnel is lower than the pool presstagee 
The bubble leaves the pore and new meniscus formation closes 
off the pore, ending this phase and returning the cycle to 
Phase i, 


Go ERP ECS OF oni. 

The introduction of oil into the pure refrigerant jee 
general, causes a decrease in the boiling coefficient as the 
Oil “concentration inereases. The effects of oil on the 
mixture are theoretically complicated as most empirical 
equations relating physical properties are derived form 
experimental data and not from theoretical equations. 
Various empirical equations are cited by Chongrungreong and 
Sauer [Ref. 5: pos.) 7084705) Yet, Jensen and Jackman 
[Ref. 16: p. 184] showed that the correlation developed by 
Chongrungreong and Sauer appeared to overpredict as oil 
concentrations increased due to poor prediction of the 
MiweurenviSscosi.cy. 

Predictive equations were developed by Jensen and 
Jackman exe ecm.o: pp. 186-187] for dense viscosities 
surface tension and specific heat for pure R-113, pure woo 


and R-113-o011]1 mixtures. These equations @2re: 
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moe 6 NC (2.4) 
Cmte 


meee Clee 7 on So (2.5) 


These equations are based on a physical model shown in 
Figure 2.4 developed by Jensen and Jackman [Ref. 16: pp. 
187-188]. The model shows a bubble growing on a heated 
surface as the refrigerant evaporates from the superheated 
liguid phase into the bubble interior. Due to the less vola- 
tile mature of the oil in the mixture, it does not evaporate 
into the bubble, but, remains at the liquid-vapor interface. 
The decrease in the rate of bubble growth is due to the 
decrease of oil diffusion into the liquid mixture and the 
decrease of refrigerant through the oil layer into the 
bubble. The decrease of the bubble growth rate decreases 


the heat-transfer rate. 


Dee LNIERNAL ENHANCEMENT 

Internal enhancement increases the heat-transfer area, 
creates additional turbulence, and produces secondary flows, 
all of which contribute to the increase in the heat 
transfer. The laminar sublayer is assumed to be the prin- 
cipal resistance to heat transfer [Ref. 17: pp. si The 
developed secondary flow'thins the sublayer and moves it 
into the stream where turbulence prevails and without an 
increase in shear; likewise, heat transfer 1S increased 
Bemenout an imerease in friction [Ref. 18: pp. SO Most 
memates show that the friction factor does increase with 
internal enhancement but the amount is strongly dependent on 
the internal enhancement geometry. The increase of both the 
friction and heat transfer is dependent on fin or rib pitch, 
groove depths, and helix angles [Refs. 17,19: pp. 6, 19-21]. 
Additionally, these studies have shown that there are 
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critical geometrical dimensions that give the maximum heat 
transfer. For example, short fins (<< tube diameter) with 
rifling will raise heat-transfer coefficients with no prob- 
lems of flow stagnation. Yet, medium fins may develop flow 
stagnation and decrease heat-transfer coefficients [Ref. 19: 
pp. 22-23). As such, the selection of internal enhancement 


will generally improve the performance of boiling surfaces. 


Z3 


IIl. EXPERIMENTAL APPARATUS 


A. OVERALL APPARATUS 

A schematic of the experimental apparatus is given in 
Figure 3.1. The photographs in Figure 3.2 show two 
different views of the apparatus. Complete details on the 
design and construction of the refrigerant and oil compo- 
nents are given by Karasbun [Ref. 20: pp. 24-32]. A discus- 
Sion of the modifications for the water heating mode is 
provided in this paper. The major components of the appa- 
ratus are: two Pyrex glass tees, an R-114 liquid reservoir, 
a water-ethylene glycol mixture sump, an R-12 refrigeration 
System, a vacuum pump, a water supply tank, three centri- 
fugal pumps, a flow meter, three heaters, a graduated oil 
cylinder, and an oil reservoir. The R-114 boiling and 
condensation occurred in the lower and upper glass’. tees, 
respectively. The R-114 vapor was condensed by the water- 
ethylene glycol mixture pumped through a copper condenser 
coil located in the upper glass tee. A 1/2-Ton R-1l2 refrig- 
eration system maintained the water-ethylene glycol mixture 
between -18 and -14 °C. | 

The major difference between the experimental apparatus 
used by Karasabun and Reilly [Refs. 20,8: pp. 24-32, 30-35], 
and the one used during the present investigation concerned 
the use of water heating instead of electric heating as a 
means to provide heat to the boiling tube. Filtered tap 
water supplied from a storage tank (15), was pumped by two 
1/2-hp motors (13 and 14) connected in Series. A three-way 
valve (V-16) was provided on the downstream side of the 


second pump to obtain two flow paths. The normal flow path 


for data ruins was Cheoupn: the metering valve (V-17), the 
flowmeter (6), the inlet mixing chamber (16), the boiling 
tube, the outlet mixing chamber (16), the chambers provided 
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with heaters (7 and 8), and returning to the storage tank 
through valve (V-20). The by-pass flow path was the same up 
to the three-way valve then the path was through heater 
chamber (8) and back to the storage tank through valve 
(V-20). Valves V-17 and V-21 were closed in order to 
isolate the boiling tube. This by-pass option was essential 
for preheating the water to higher temperatures, to achieve 
higher heat flux values in the boiling tube compared to 
values obtained from room-temperature water. A 500-W and two 
lO00-W heaters were used to maintain a steady water inlet 
temperature and to preheat the water. Flexible pressure hose 
was used for water piping throughout the apparatus. The 
mixing chambers and sections of hose connecting the mixing 
chambers to the boiling tube were double insulated to ensure 
SaculbraLe temperature measurements. 

Copper-constantan thermocouples were used to measure 
thermal emf's for R-114 Fanta @and vapor, water inlet and 
outlet, and water-ethylene glycol mixture. A series- 
connected thermopile consisting of 10 junctions on either 
end was used to meaSure the temperature drop of the water 
from the inlet to the outlet of the boiling tube. This ther- 
mopile was calibrated against two quartz thermometer probes 


and agreement was found to be better than + 0.02 K. 


B. BOILING TUBE CONSTRUCTION 
Figure 3.3 Shows schematics of the smooth tube, and 
two enhanced tubes used. The tubes tested were: 

l. a smooth copper tube; 

2. a _.90:10 copper-nickel corrugated tube (commercially 
referred to as Korodense tube with an external, 
Sintered porous coating (i.e. High Flux); 
an internally and externally enhanced tube (alloy 
eZ 200). Beaucee Dyewolvermumne Division of U.Q0.P. 
(commerically referred to as Turbo-B); 


4. a porous-coated (High Flux) Korodense tube with the 
porous coating machined off; 


on eee tube with the external enhancement machined 
off. 
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These tubes will be referred to as High Flux, TurboeB ieee 
fied High Flux, and modified Turbo-B tubes, respectively, 
with the exception of tubes 2 and 4. - Tubes 2 and 4 will be 
referred to as Korodense tubes with respect to discussions 
of only the internal enhancement. All of the tubes had an 
active (i.e. , heated) section of 304.8 mm in length. The 
smooth, High Flux, and modified High Flux tubes had an outer 


diameter of 15.9 mm and an inside diameter of 12.7 mm. in 


‘the case of the High Flux tube, these diameters represent 
the values in the smooth portion of the tube (i.e., before 
performing the corrugation process). The Turbo-B tube had 


an external structured surface and a multiple-start helical 
ridged inner surface. This tube had an outer diameter of 
16.9 mm at the base of the external structured surface and a 
minimum diameter of 14.5 mm at the internal ridge tips. The 
modified Turbo-B tube had identical diameters as the Turbo-B 
tube. The 63.5-mm-long inactive sections on either end was 
insulated by Teflon sleeves of 1.6 mm wall thickness located 
inside the boiling tubes, as shown in Fieure foe The 
requirement for using these two modified tubes is explained 
further in Chapter IV (EXPERIMENTAL PROCEDURES). 


C. COMPUTER-CONTROLLED DATA ACQUISITION AND REDUCTION 

Hewlett-Packard equipment was used for data acquisi- 
talon. control, and analysis, as well as for sStoimoeeon 
data. An HP-3497A acquisition unit was used to read the 
thermocouple and thermopile outputs. Table 1 gives the 
channel designations used. 

Information entered by keyboard to an HP-9826A computer 
unit prompted and controlled the data acquis vem unate 
Data were analyzed and stored with the same computer unit. A 
Step-by-step description of the data-reduction procedure is 
given in Appendix A, along with a printout of the prograin 
(DRP4) used. 
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TABLE 1 
CHANNEL DESIGNATIONS 


Channel Function 
~ R- iS a thermal emf’s T(1) 


2 R- 114 ae thermal emf T(3) 
3 water-ethylene glycol mixture 
thermal emf T(4 
4 water inlet thermal emf T a 
5 water outlet thermal emf ) 
20 thermal emf difference across 


water inlet and outlet T(7) 
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IV. EXPERIMENTAL PROCEDURES 


A.  FLOWMETER CALIBRATION 

This section gives a brief description of the flowmeter 
Calibratsom: a more detailed explaination is given in 
Appendix B. A Fischer Porter flowmeter was used to indicate 
percent water flow through the boiling tube. Prior ee 
conducting any boiling runs, it was necessary to develop a 
correlation relating the flow percent to the mass flow rate 
and water velocity. Additionally, any temperature effect on 
viscosity would be incorporated into this correlation as a 
COrreckmon Lactor. Calibration runs were conducted over a 
temperature range of 19 “C to 38 °C and varying flowmeter 
Settings. The correlation (equation (B.1)) was determined 
from these results. The correction factor for the viscosity 
temperature dependency was determined to be 1.0 with an 
accuracy of + 0.02 as explained in Appendix B. In fact ieee 
was not possible to find a systematic trend for this correc- 
tion factor with the water inlet temperature. The stated 
accuracy mainly consists of the uncertainties involved with 


the visual reading of the flowmeter settings. 


B. BOILING TUBE AND APPARATUS PREPARATION 

The external surface of the boiling tube was’ cleaned 
with a Nitol (2% nitric acid and 98% ethyl alcohol) sole 
and both the external and internal surfaces were rinsed with 
acetone and were air dried prior to installation in the 
apparatus. A vacuum test at an absolute pressure of about 50 
mm Hg was conducted on the refrigerant side of the apparatus 
after the boiling tube was installed. If no signs of leaks 
were evident after two hours, the system was pressurized to 
a gage pressure of about 250 mm Hg with R-114 vapor and the 


apparatus was checked for leaks with an automatic Halogen 
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beak Detector. After fixing any leaks that may have been 
present, R-114 liquid was transferred from the reservoir to 
a pre-determined level in the boiling tee. This level gave a 
free surface 40 mm above the centerline of the boiling tube. 
At the set level, the mass of the R-114 liquid was computed 
to be 2.48 kg. A desired saturation temperature of 4 °C and 
the resultant system pressure were obtained by turning on 
the R-12 refrigeration system and varying the flow of cold 


(approximately -16 cc) water-ethylene glycol mixture through 


the condenser coil with valve V-9. At this saturation 
temperature and pressure, each boiling tube developed 
nucleation sites. The saturation temperature of 4 "Cawas 


maintained for at least 30 minutes prior to initiating 
nucleate boiling on the complete active length of the 
boiling tube. 

The initiation of nucleate boiling could be done with 
cold water at 19 °C or with warm water at 25 °C through 
valve V-17 or with steam from a steam generation system 
through valve V-22. Early experimental data showed slight 
differences (up to 10 percent) in the computed outside heat- 
transfer coefficient depending on the highest heat flux at 
which the boiling was initiated on each boiling tube. As 
will be discussed in Chapter V (RESULTS AND DISCUSSION), the 
steam intiation method provided the highest initial wall 
superheat and ensured the beSt possible repeatability of the 
data. This intiation method was used for all of the runs 
unless otherwise specified. 

Steam initiation was conducted by introducing steam 
through valve V-22 and discharging steam drain valve V-18 
for about one minute. After this period, the three-way valve 
V-16 was shifted from a neutral position to the normal flow 
position. Valves V-22 and V-18 were closed and valve V-21 
was opened simultaneously as the three-way valve was shifted 


to the normal position. Valve V-17 was opened wide prior to 
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Shifting of the three-way valve. As these valves were shift- 
ed, warm water at 25 °C was pumped through the boiling tube 
at the maximum water velocity, which in turn gave the maxi- 


mum heat flux for the given water inlet temperature. 


C. BOILING DATA RUNS 

Once the nucleate boiling initiation was completed, the 
saturation temperature was raised to 13.8 °C by adjusting 
valve V-9. This valve was also used to maintain his Ssatuiee 
tion temperature within +0.05 K throughout a boiling run. 
The water inlet temperature was maintained to withiny2070Gem 
by using the heaters or introducing filtered cold water, 
respectively. The combination of heaters and cold water was 
used to maintain the desired water inlet temperature during 
each run. Two types of boiling runs were conducted at each 
oil mass percent. 

The first type of run was conducted after holding the 
desired saturation temperature, water inlet temperature and 
initial flowmeter setting at steady-state conditions for 
about 10 minutes. After this period, the flowmeter setting 
was decreased with two data sets being taken at each flow 


setting. The two data sets were to show repeatability at 


each setting. The period between pairs of data sets was 
about 5 to 8 minutes, with steady state for each flowmeter 
setting maintained for at least 2 minutes. The second type 


of boiling run consisted of maintaining the previously 
mentioned initial steady-state conditions at a selected high 
water inlet temperature. The run was commenced by intro- 
ducing cold water into the storage tank which lowered the 
water inlet temperature. While holding the water velocity 
constant, six to seven data points were taken, on a one time 
pass, at about every 0.3 K inlet temperature decrease. When 
the inlet temperature reached a selected low temperature, 


the steady-state conditions of the low inlet temperature, 
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water velocity, and saturation temperature of 13.8 °C were 
maintained for about 10 minutes. Then, the water inlet 
temperature was increased using the heaters, while holding 
the water velocity constant. Again, data points were taken 
on a one time pass, in about 0.3 K intervals, ieee 
inlet temperature reached the previously used high inlet 
memperature. 

Two complete boiling runs of the first type mentioned 
above were performed at 0% and 24 oil concentrations’ to 
further demonstrate repeatability. Table 2 shows the percent 
difference in wall superheat (i.e., AT) and heat flux for 
each boiling tube at these two oil percents. The higher 
percentages for the High Flux and Turbo-B tubes for the wall 
Superheats are Simply because of the low wall superheat 
values of these two tubes when compared to the values for 
the smooth tube. At these low superheat values, a small 
difference in the wall superheat leads to a larger 
percentage difference than that which occurs at higher wall 
Superheats. The variation of the water inlet temperature was 
not considered a factor in this percent difference as the 
variation was less than + 0.02 K per minute during these 


constant-inlet-temperature, decreasing-water-velocity runs. 


TABLE 2 
BOILING RUN REPEATABILITY 


Tube al 7Z ok q 
Smooth oy 70) 
High Flux 0 £309 £194 
Turbo-B 0 +307 a 47 
Smooth 2 +107 DE 
High Flux 2 +102 z ay) 
Turbo-B +307 + 2% 


Note: initiations performed with steam 
Additionally, data runs were conducted to determine the 


temperature increase across the boiling tube due to the 


pressure drop between the locations of the thermopile 


oS 


probes. The data collected showed that this difference was 
less than the +0.02 K accuracy of the thermopi lean. 
considered negligible. 

Following the boiling runs in pure R-114, oil was 
introduced to the boiling tee through valve V-l from the 
graduated cylinder (5) Shown in Figure 3.1. After the 
required volume of oil for a given oil mass’ percent was 
added, water at the maximum velocity was pumped through the 
boiling tube to promote vigorous’ boiling. This boi Fagg 
ensured good mixing of the R-114 and oil mixture. Either 
another boiling run was conducted or the system was shut 


down to prepare for another steam initiation. 


D. DATA ACQUISITION AND REDUCTION 

The automatic data acquisition system was prompted to 
record the required thermal emf's, by keyboard inputs to the 
computer. The data were immediately processed and printed 
out on a hard-copy printer. A step-by-step procedure of data 
processing and reprocessing 1S given in Appendix A. Also, 
a printout of the data reduction program, DRP4, is included 
in this appendix. The heat flux, overall heat-transfer coef- 
ficient, outside boiling heat-transfer coefficient, and wall 
Superheat values were computed based on the outside area 
expressed by the diameter if the enhancement of the outer 
surface was removed. Further, in order to account for the 
heat conduction through the inactive end sections, a correc- 
tion was made for heat flux using an iterative computation 
procedure based on natural convection as discussed by 
Karasbun [Ref. 20: pp. 54-56]. 

When a boiling run was completed, the data were repro- 
cessed by computer to obtain the inside and outside heat- 
transfer COCEFICLTENES Ons further calculat tons: A 
Sieder-Tate-type constant (C,;) for the inside coefficient 
for each of the externally smooth tubes was calculated by 
uSing a modified Wilson plot | Ref. 22) . For the smooth tube, 
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modified High Flux, and modified Turbo-B tubes, the outside 
EQeiricient was calculated simultaneously using the cor- 
relation (equation: (e355 and (C.4)) developed by Rohsenow 
beet. 23:  p. 969], and the modified Wilson plot. The 
Rohsenow correlation is based on externally smooth tubes, 
thus the requirement for the modified High Flux and Turbo-B 
tubes. A further explaination of the modified Wilson plot 
is given in Appendix C. ieee aluecetound tom the sexter- 
nally smooth tubes were used for the corresponding exter- 
nally enhanced tubes. Knowing the inside coefficient, the 
outside coefficient for -the boiling tubes were computed by 
subtracting the inside and wall resistances from the meas- 


ured overall thermal resistance (equation: (C.7)). 
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V. RESULTS AND DISCUSSION 


A. INSIDE HEAT-TRANSFER COERFREIGIENDT 

Results were obtained for the smooth, High Plo and 
Turbo-B boiling tubes at a boiling temperature of 13.8 ac 
and at oil concentrations of 0, 1, 2, and 6 mass percent. 
The heat flux versus wall Superheat for these three tubes is 
Shown in Figure 5.1. All three tubes are in the nucleate 
boiling range with decreasing heat flux. It can be seen that 
the High Flux and Turbo-B tubes outperform the smooth tube 
throughout the tested heat flux range. Additionally, the 
Rohsenow correlation for a smooth tube at a boiling tempera- 
ture of 13.8 °C with the experimentally determined (from the 
modified Wilson plot) Cgy coefficient of 0.0060 is shown in 
this figure. 

The Rohsenow correlation, as used in the modified Wilson 
plot, has an exponent of r-= 1/3. Yet, this correlation is 
very sensitive to the r value. A change in this exponent of 
three percent yielded a change in the Cos only of 0.016 
percent but gave a change Of Ly bercent in the 
Sieder-Tate-type constant inferred from the modified Wilson 
plot. For the purpose of data analysis, the value or r = 1/3 
was used, while knowing that an additional uncertainty was 
being introduced by this exponent. 

For each externally smooth tube, the Sieder-Tate-type 
constant was computed using the modified Wilson plot based 
on 5 to 10 different runs in pure R-114. The computed 
values for the smooth, Korodense, ‘and Turbo-B_ tubes were 
0.036 £0.002, 0.066 + 0.005, and 0.077 +t 0.003, resPee. 
tively. For the smooth tube, the experimentally found 
Sieder-Tate-type constant is 33 percent greater known than 
the well-known value of 0.027 for long internally smooth 


tubes with fully developed _ flow. This larger constant 
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appears to be the result of the entrance effects of the 
shorter experimental tube and the uncertainty introduced by 
the exponent r of the Rohsenow correlation as discussed 
above. The minimum entrance length for a fully developed 
pipe flow is given, as a rule of thumb, by Incropera and 
DeWitt [Ref. 24: p. 406] as 


L 
AS x 60 (Ss. 
D 

for Reynolds numbers greater than or equal to 10000. Also, 


as the pipe roughness increases the minimum entrance length 
decreases. : 

Withers has developed Stanton mumber correlations for 
tubes with single-helix and multiple helix Tidgume 
[Refs. 25,26: pp. 52-56,44]. These correlations are listed 
below.” ~ For single-neli<sc id sine. 








(5.29 
(5.39) 
where r and m are determined by tube dimensions, and 
ee {2.540 P(g) + 3.75 (5.4) 
For muleiple-helim cideine. 
—_—_ 7 - 
Se 5 Se ee : (5.53 
: f a 
8 ( Ref = + ¥Y 
og l 
| 7 pa a a (5.6) 
oo +Oxn | T Re / | 


rand m again are functions GE tube Gimemomonc. 
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A comparison of Stanton numbers was made for the modified 
Korodense and modified Turbo-B tubes using the following 


equations for these two tubes: 





ns De 
st = o—_—_——- (525) 
eee Pr 
where 
k 0.3 Oss. Gee Oe: 1 4 
4 = Ge eee ee ae 
a i o0; 7" ( - (33 0)) 


Mee @eeStiits for these two modified tubes were compared to a 
copper Korodense ( er, Single-helix ridging) and a 
Turbo-Chill (i.e., multiple-helix ridging) tubes, respec- 
tively. The latter two tubes tested by Withers were similar 
in radial dimensions but were 1.5 meter long. The predicted 
St numbers using equations (5.1) -(5.8) for the Withers’ 
Korodense and Turbo-Chill tubes were 0.00177 and 0.00273, 
respectively. The experimental St numbers using equations (5.9) 
and (5.10) for the modified Korodense and modified Turbo-3 
tubes were 0.00242 and 0.00286, respectively. The majority 
of the difference between the predicted and experimental St 
numbers may be attributable to the entrance effects and the 


Mncertainity introduced by the Rohsenow correlation. 


Bee ULGHT OFF EFFECTS 

AS referred to in Chapter IV, Figure 5.2 shows the light 
@eewetiects for the smooth, High Flux (i.e., porous-coated 
Korodense), and Turbo-B tubes. These effects are most prob- 


ably due to the requirement of greater heat flux to initiate 
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Eneminne Leation Sites than the heat flux required to maintain 
the sites once they are activated. As the boiling tube sits 
in the R-114 liquid, prior to initiation, reentrant cavities 
om the mucilecation sites flood with lagquid. The initial 
formation of the vapor bubbles and thus the initiation of 
the nucleation site requires the input of a greater heat 
flux or higher wall superheat. Figure 5.3 shows the boiling 
tube imitiation with steam for the smooth, High Flux, and 
Turbo-B tubes [Refs. 21,7: p. 478, 62]. As shown in Figure 
ao, ticmetonerii tube 15 more sensitive to laght off 
effects. This sensitivity is mostly due to the lower wall 
Superheats at which nucleation occurs. A change of .5 of a 
degree iS more pronounced at these lower wall superheats 
than a similar change at the wall superheats for the Turbo-B 


and smooth tubes. 


e. SMOGHE TUBE 

Figure 5.4 shows the performance of the smooth tube at 
wee Ll, 2a ana 6 mass pereent oi! concentration at 13.8 Ge 
Reilly's [Ref. 8: p. 65] data for the electrically heated 
smooth tube at a boiling temperature of -2.2 °C in pure 
R-114 liquid are also included for comparison. The differ- 
ences between the experimental results shown and Reilly's 
data are probably due to physical characteristics introduced 
in the electrically heated tube when soldering the copper 
sleeves (inside of which the electric heaters were fitted) 
on the interior of the tube. As discussed by Reilly [Ref. 8: 
pee /-Giliewetnmewcomtact resistance at the interface between 
the sleeve and the inner surface of the tube was minimized 
by tinning. However, the tinning process may not have been 
100 percent successful, thereby introducing unacceptable 
uncertainties into the experimental meaSurements. Further 
Reilly reported considerable variations (up to 3 K at a heat 
flux of 98 kW/m2) in the measured wall temperatures. Reilly 
attributed this observation to an axially non-uniform heat 


flux generated by the cartridge heater. 
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As shown in Figure 5.4, the wall superheat increased for 
1 percent oil concentration and decreased slightly for 2 and 
6 percent oil. This decrease in the wall superheat for the 
higher oil percentage is possibly due to the enhanced bubble 
formation experienced by the tube because of the foaming 
action of the R-114 and oil mixture decreasing resistance to 
heat transfer as explained by Chongrungreong [Ref. 5: pp. 
701], Nobukatsu [Ref. 27: p. 60], and Chaddock | Reiy 21505 
474]. 


D. HIGH FLUX TUBE 

The performance of the High Flux tube in terms of heat 
flux versus wall Superheat 1S given in Figure Syoe 
Additionally, this figure shows Reilly's data [Ref. 8: pp. 
70], for the electrically heated High Flux tube in pure 
R-114 at a 6.7 °C boiling temperature. The presence of oil 
increased the wall superheat, at a heat flux of 40 kW/m*, by 
a factor of 1.22 for 1 percent oil and a factor of 2.12 tae 
6 percent oil. The number of nucleation Sites reduced 
considerably and the foaming action increased significantly 
as the oil mass percent increased. TheSe actions are shown 
in Figure 5.6. With 6 percent o1!] comecentratemern, it was 
visually observed that the nucleation sites were active only 
in the corrugations, as shown in)’ Figure 95,6 (ce - This 
congregation of sites is due, in part, to the sparser porous 
coating on the high points of the tube and the decreased 
wall thickness in the corrugations. As a result of the 
corrugation process during tube manufacture, the wall thick- 
ness at the corrugation is smaller than the rest of the 
tube, which implies lower wall resistance. During the 
Sintering process, the coating appears to have concentrated 
in these corrugations giving a thicker coating of coppem 
particles and increased number of nucleation sites compared 


to the rest of the tube surface. 
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The increase of the wall superheat, when oil 1s present, 
is possibly the result of increased surface tension and 
viscosity. The increasing surface tension and viscosity lead 
to the formation of an oil-rich layer next to the boiling 
surface. This oil-rich layer has an insulating and inhib- 
iting of bubble growth effect within the cavities. The 
higher liquid-vapor surface tension and surface tension in 
the oil-rich layer requires more energy for the formation 
and growth of the bubble. Additionally, diffusion processes 
affect the bubble formation. Oil concentrations near the 
boiling surface are high compared to the concentration in 
the bulk liquid and diffusion occurs from the cavity into 
the bulk liquid. This diffusion continues until an equilib- 
rium point is reached. If oil-rich layers are developing on 
the boiling surfaces, the diffusion from the cavities to the 
bulk liquid slows down aiding the further build-up of the 
oil-rich layer. The vapor bubble forming in the cavity also 
has an oil vapor content which is undergoing its own diffu- 
Slon process. As the refrigerant evaporates through the 
Oll-rich layer into the interior of the bubble, it must 
overcome large diffusion resistances set up by the oil vapor 
within the bubble. The increased surface tension, viscosity, 
diffusion resistance and the decreasing bubble formation 
rate contribute to the decrease in the heat-transfer coeffi- 
clent and the increase in the wall superheat at a given heat 
Mme kets. 5,21,238929730: p. 702, 477-478,59,372, 82] 

Also, boiling runs were conducted holding the water vel- 
ocity constant and allowing the water inlet temperature to 
decrease and then to increase. The decreasing and increasing 
heat flux with the respective change in the water inlet 
temperature versus wall superheat is shown in Figure 5./7. 
The motion of the hysteresis (i.e., clockwise or counter- 
clockwise) depended on the temperature starting point. If 


the temperature was decreased and then increased, the motion 
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was clockwise. When the temperature was increased and then 
decreased, the motion was counterclockwise. This directional 
change is related to the amount of heat flux required to 
maintain active nucleation sites. In the first case of tem- 
perature decreasing then increasing, less heat flux is need- 
ed to maintain the same amount of sites, but, at the same 
time, heat flux decrease as the water heat capacaity 1s 
decreasing. In the second case, with the water inlet temper- 
ature increasing then decreasing, the water heat capacity 
increases so the heat flux increases raising the wall 


Superheat and the heat-transfer coefficient. 


eae 6LLEURBO-B TUBE 

Figure 5.8 shows a similar relationship as Figure 5.5 
but for the Turbo-B tube. The increase in wall superheat, 
for a heat flux of 40 kW/m*, ranges from a factor of 1.04 
mor - percent o1l to a factor of 1.37 for 6 percent oil. The 
Turbo-B tube remained active over the full active length of 
the tube, as shown in Figure 5.9. The apparent decrease of 
nucleation sites, as shown by the High Flux tube, did not 
occur with the Turbo-B tube. As discussed pertaining to the 
High Flux tube, the increase in the wall superheat with 
increasing oil concentration is probably due to the 
increased surface tension and viscosity of the mixture 
Within the cavities. Also, Boiling runs for constant water 
velocity and varying water temperature were performed on the 
Turbo-B tube and are also shown in Figure 5.10. The same 
hysteresis motions observed with the High Flux tube are 


observed with the Turbo-B tube. 


F. COMPARISON OF BOILING HEAT-TRANSFER COEFFICIENTS 

The ratios shown in Figure 5.11 are h/h for the 
Turbo-B and High Flux tubes and h/h, for the smooth tube. 
All of the ratios were taken at a heat flux of 40 kW/m?. The 


High Flux tube shows the most significant effect due to oil 


Dik 


“IQNL FT-oOqiny 107 
qeoyszedns TTCM YIN XNPF Away JO VUoOTIvTIvA “Qs DINIT I 


On @O'h-= 1) 





rote x 9 ,OT X S‘T 
ot X € 
pueq Lazyuypearasun + O41 
Oe, 3 
to LY 2% —— — Lo 7 
Te 2 ee 
tro a 
to LOT) / vy 
AY 
Wa 
AG n cl! 

{Nn 
x. 
~ 
4 
nN 
VV 





—— m0 X 6 


6) Crd 4 


Cowra 


O 


” 


A 


gs 


_— ea 
oe ee 


he 


CAV 


Cw 








“Ole xa 9 


agny 
rajyey Surpsuvyy Yat 


Piuyerrssrep Jwry 


AOJ VIN ArAIdWIL, JOLTUT 


XNPA WY Jo 


O) (2 i= ab 


—= 


NOt rae A oe) tS 


(Buisearsag way 
Mursrzazouyg dury 


E06 9 oa= 
110 7% 
TtOe: 


Tto % 


oT) 


Sq 


‘11 4&5 tty) 





,Ot X 8 


eB TOF YuTIIIg ssep [TO YAM OTArY 
Taysuerp_-jeoyy JUTTIOG Jo wuoTArTArA 


yIoows 


q-oqinr © 


xOTs UFTH 


V 


"su/M OF FO XnNTA Wey 


qRUPOUBY 110 
v 


"IT'S aandstTy 


48 


ca 


GT 


Et 


YUyzooueYy/y 


55 


concentration which is due to the reasons previously 
discussed. The h/h, ratio for the High Flux tube decreages 
by 65 percent. For the Turbo-B tube the rate of degradation 
of the boiling heat-transfer coefficient is less pronounced. 
The Turbo-B tube shows a decrease of 24 percent at 6 mass 
percent oil. The larger cavity size of the Turbo-B tube, 
compared to the High Flux tube, results in smaller oil 
concentrations in the oil rich layer. This observation may 


explain the smaller degradation experienced by Turbo-B tube. 


G. OVERALL HEAT-TRANSFER CHARACTERISTICS 

The overall heat-transfer coefficient versus decreasing 
water velocity for oil mass percents of 0, 2, and 6 for each 
boiling tube is shown in Figure 5.12. For pure R-114 with a 
water velocity of 2m/s, the Turbo-B’ tube outperforms the 
High Flux tube by a factor of 1.13 and the smooth tube by a 
factor of 4. While the High Flux tube outperformed imme 
Turbo-B tube based on the outside heat-transfer coefficient, 
the reverse is true for the overall heat-transfer coeffi- 
CLent. TaissGelartilonship — am. U, is due in part to’ the 
increased internal enhancement of the Turbo-B tube. 
Comparison of the Sieder-Tate-type constants for the Turbo-B 
and High Flux tubes (see p. 34) shows that the Turbo-B tube 
has a 1/7 percent greater inside coefficient than the High 
Flux tube. Both of the enhanced tubes dramatically outper- 
formed the smooth tube. It should be noted that the increase 
in pressure drop due to internal enhancement (compared to a 
smooth-interior case) was not considered here. The pressure 
drop would be an important factor in finally selecting the 


most economical tube type. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


Considering the data gathered in this investigation for 


boiling of R-114 and oil mixtures at 13.8 °C boiling tempers 


ature, 
ime 


the following conclusions are reached: 


Based on the overall heat-transfer coefficient, the 
Turbo-B and High Flux tubes outperformed the smooth 
tube by a factor preater than 3, and the Turbo-B tube 
outperformed the High Flux tube by a factor of : 


At a heat flux of 40 kW/m* with zero mass percent oil 
the outside coefficients of the High Flux and Turbes 
tubes were factors of 14.6 and 6.4 times, respec- 
tively, compared to the smooth tube. These factors 
decreased to 7.0 and 4.9 for the High Flux and Turbo-B 
tubes, respectively, with 6 mass percent oil. 


The Turbo-B and High _Flux tubes showed Sieder-Tate- 
type constants of 0.0/7 and 0.066, Cee oe aes com- 
pared to the Sieder-Tate constant of 0.02 for long 
smooth tubes. As noted in_the discussion, these values 
may be up to 10 pee lower for internally enhanced 
long tubes than the values for the short tubes. 


While the High Flux tube outperformed the _Turbo-B 
tube, based on the outside heat-transfer coefficient, 
oe High Flux tube is more susceptible to the presence 
Om onl. 


B. RECOMMENDATIONS 


Based on the results obtained from this investigation, 


the following recommendations are made: 


IS 


The investigation of water-heated tubes should be 
Ones ace to use long tubes. This would decrease the 
entrance effects that short. tubes experience and it 
would come closSer to duplicating actual operating 
conditions. Also, the pressure drop effect prev ona 
mentioned in Chapter V (RESULTS AND DISCUSSION could 
be investigated in the analysis of the long tubes. 


Studies of the boilin heat-transfer coefficients 
Should be expanded to tube bundles. The interaction of 
tubes within the tube bundle is a major factor in the 
analysis of heat-transfer data. The 1a heat- 
transfer coetiicient ny, be strongly influenced by any 
foaming action within the bundle. 


Further investigations should be conducted on internal 
tube enhancements due to the major thermal resistance 
being the internal resistance. 


Modifications to the condensing side of the present 
apparatus would allow for studies of short water- 
heated tubes at lower boiling temperatures. The major 
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modification would involve only the increase of the 
condenser drain piping to handle increases in the flow 
rate of the refrigerant liquid. 


a 


APPENDIX A 
DATA REDUCTION PROGRAM 


The data reduction program, DRP4, used for this investi- 
gation is listed below. A brief description of the major 
Sections iS given in Table 3, which is followed by an indi- 
vidual line listing of DRP4. 
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Line Number 


10-30 
S0-3070 


3700-4165 


4170-4220 
a2 253935 


4940-6175 
6180-6410 
6415-6490 
6495-7280 


Loo 7 370 
7395-9640 
9645-9725 


9730-97 52 


77 OS OZ 00 


i [OAS = 108 1o0: 
10365-10475 


TABLE 3 
DRP4 MAJOR SECTIONS 


Description 


Selection options (i.e. process 

Gatiaw plo, etc. 

Sub Program Main 

-boiling tube dimensions 

-select desired ae 
conditions (i.e. initiation 
Node sulowerate., Fsat, ete: 
-convert emf s to temperatures 
[~COMpPULe CONEACE TeSi1Stanee 
(electrically heated mode only) 
,cOMpUubre VaGlOUs Walter properties 
-various heat transfer 


ealemlations(@, EMED, U etc. ) 
-compute various Freon 114 
PrEOperetes 


-compute natural convective heat- 
EGamoter cOetitecLtent for 
unenhanced ends(iterative 
procedure ) 

ecoupute heat loss for unenhanced 
ends 

-compute actual heat flux and 
boiling coefficient 

-record and store data 

RUMGEDOnNS EOteecurve fits of 

various Freon 114 properties 

Rimetion) fOrepolynomial curve fit 

suo Prorram Fowy ; 

-calculate nth order polynomial 

sub Program Plin 

-plot linear eee 

Sub Program Stats 

-calculate standard deviations 

Sub Program Coeff. 

-plot cross-plot files 

Sub Program Wilson 

-calculate modified Wilson 
plot values (See Appx. C 

Pumetions for curve fits of 

various water properties 

sub Program Plot 

-plot non-linear graphs 

sub Program Symb 

-print various symbols on 
plotted graphs 

Sub Program Purge. 

-purge unwanted files 

Supe eeoecrcam Ldcn 

-calculate water temperature 
rise corrections due to 

ressure drop between 
hermopile locations 

Sub Erorram U plot 

Sjeheakgue files 

sub Program Select 

-select various sub 
program options 


ok 
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oO PRINTe? IS 
4) CALL Sereer 
a5 INPUT “NANT TO SELECT ANOTHER OPTION (1=%.9=N)7", Tsai 
$0 brs teei=! THENSGOIO 40. 
Ss Beas 
50 BES 


55 OR eRe 

74 PRINT “DATA COLLECTION/REPROCESSING COMPLETED” 
‘5 END 

30 Suse ain 

35 COM tas7 male 

39 COM /Cc/ CC7 ye teal 

a5 COM SW? D2. Di Dew ce Vitae 

190 DIM EmfCl2).7TC12) .dlata) .D2atS) 01at 3) Deal 3), Lata) buat 3) e026 eee 
(AS cap S| 

mS VATA 2.7008609! .25727.34282 767545762 )5.7 0025595751 


150 DATA =9324/7485589,5.37648E+11,-2,65 1 392E~+13 .3.354078E+14 
115 RESD Cee) 7 

120 ' DATA “Smooth”. Hian Flux”. “Thermoaexel-©”.”°Thermoexel-—HE - 
> DATA Smooth.Hign Flux.turoo-8.Hign Flux Mod.Turoo-8 Mod 
je READ TnS¢s) 

TO PRgieR 6S - 708 

140 Beer 


fas IF Idp=4 THEN PRINTER [S ! 

1$0 TF fap=4 THEN GOTO 1280 

P5 IMPUT “ENTER MONTH. DATE AND TIME (MM: DOD:H4H:4N:SS)",Dates 
160 OUR 709. Dases 

165 QUTIPUT 709:°°TD” 

170 ENG 022Dates 

(75 PRT 

180 PRINT ” Month. datceranaeclmer<: (Dares 
185 PRINT 

ea Pell JOIlG “10k. “Ie Peceranm pane Gaon s 
ees Bees 

CNY LNPYUT “EMTER Die MUMBER On 

205 PROAE USING. “16x% Cisx ontmoee = ee in 


a 


20) SE ar 

oe Lue Eee SUPUY AO0e Uo) eee 

20 Bee 

oe Peet» Sele HEATING MODE (O=ELEC;, t=NATERD”. Ihm 


med BES. 
Eo tpt “eslteR “RESMOCOUPLE TYRE CieNp 6; ="LD so eal 
e 


241) Lr im=)) ie N 

os le! 

oe INPUT “GTYE A NAME FOR THE RAR DATA FILE” .D2 rites 
255 PRINT USING "16X.""New file name: "".14A":D2_ es 


259 Sea 
25 CREATE 3DAT D2 1 
nay ASSIGN @wFiiee TO D 


-— 


foW 


aol 


uy t(D 


iT ts 


29 ' DUMMY FILE UNTIL Nrun <NOWM 
395 ot fy la s="DUMMY" 
230 CREATE 3DAT Di _ +filesS.Size! 

mF 


225 Aso 1.01 fe) 10 De es 
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340 
245 
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3oD 
360 
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330 
335 
400 
405 
4410 
At 
420 


425 


430 
43 
440 
4a5 
450 ! 
455 
460 
455 


it 


WIA NAVI 
Up Geo) i Oi 
mu 


sy ee en 


NUMe aoe we mets 17E Co (Sper 3uL 7). 


re) 


ePuy 
[ee iadte= 
Late! =) 
Pate. =] 
Beef) JothG “'6A, «lo detective 'Cs exist” 
END TF 

ieee c= | 
Baoe 
INPUT “ENTER DEFECTIVE TC LOCATION’ .Ldtc! 

BeepmmuotNG “lok, JC is eereetive at iocation °".D”:Ldtc! 

aero 2 =() 

END IF 

Marldte2s2 HEN 

Bee - 

elie so Mens Deeacr VE TC SiGsT ees iedtc! .Ldtc2 

eee Sieve C aremdeqoetive at locations ~'.0.4X.0°:Ldte!l .Ldtc2 
ND i 

TF Tdtc>2 THEN 

Bear 

Bom iieRo lS | 

BEEP 

PRINT “INVALID ENTRY" 

PRINTER IS 701 


iHEN 


Gata si0 

ee SIF 

END IF 

CMe) BF rle::Lote! .Ldtc2 
Im=' option 

cue 

cE: 


eT “GIVE THE NAME OF TRE EXISTING DATA FILE”. D2_Filed 
PRINT USING “16X.""OQid file name: "".14A":D2 files 

SOoIGN 3File2 10 D2_-:1les 

Ses Filed: Nrun 
ele x “fF :iez:Do0ldas 
eet USING lek, this Qatawvcecrcacem on. «!4A 2De1d5 
a ecer Plier: :! dtc? ,Latc. 

pees tc! Oo) OR Latcl># THEN 

Bees iON 1oA, ihernocoupies were derective at locations: ”.2(3D.4x%)"; 
Teaco. 

Spe Tr 

see We ple2elic 

ED l= 

te Im=) AND fhm=! THEN 1595 


BEL? 

Peete aetlil 10 CREATE 4 PLOT FILE? (O=N.!=Yo"  Iplor 
fF tplot=l THES 

BE=2 


eee oe eMeee WR PLOT FILE’ .P rries 
CREATE SDAI P F¥:ie5.4 
ASSIGN w¢Plot ee 


“ND IF 

Tr Ihm=! TREN 

BEES 

ENP Se ieeneeior tice COs 5 = )"" .cuct 
Bee tint=" TREN 

BEe> 
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‘ee 
390 
305 


339 


B15 
330 
325 
opel 
295 
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7 oo ow. 


LPH “SNES go. = k= “AME”. Pp oes 
Cheat] 2p =e oe 
aco Gul GS Yt eS 
eHD [fF 
Sees 
INPUT “SANT TQ CREATE Re FILE? (GaN, t=)" Tee 
IF tre=1 THEN 
RED 
INPUT “ENTER Re FILE NAME” .Rerirlesd 
CREATE BDAT Reriies.10 
ASSIGN w%Rerile TO Reriles 
END IF 
ENO eLe 
PRINTER IS 1 
IF Im=0 THEN 
BES? 
PRINT USING “SX.° "Select tube numoer ” 
ita Thm=9 THEN 
PRINT USING, “Sk. "0 “Smooth 4imen mer 
PRINT USING “SX.°°7 Smooth 4 inen Cu (Press/sciae 
PRINT. USiibe 64. 2° “Sort sou tomremines eo ee 
PRINT USING)” 5k.°°3 -Sott Solder <4 unch BIGH Fee 
PRINT USING "6X."°"4 Aiteland Hard 3 inch’”” 
PRINTSUSING “6X6 "SS: CHIGHIr EU aeernen 
PRINT USING "6A. "5 GEHA-—=K 13 Fins/in. — 
PRINT USING “BX.°"7  GEWA=K 25 Fins/in 
PRinkt UsinG.= 5A\ fe Genhe |) Too inev on: 
PRINT USING "6X.°"°S GEWA-T 26 Fins/:in’”” 
Else 
PRINT USING "“SXK.°"9 Smooth tube’ 
PRINT USTHGs Gh. 1” “Higmer lice 
PRI USiNGtcabA, 2 Mitoo-S oe 
PRINT USING “6X. "3S High Fluxeroa ~~ 
PRINT USING “S4%.°°4 “Turge=8 Hod” 
END IF 
PUP bee 
BUTE wear lelelts 
See 
PRINTER SIS 7 00 
TF fetcig THEN PRIMT USING “164A. | ube *hinber: Tee lata 
TP fetoo THEN PRinT USING "T6X, ibe Number : A) Deel ec 
iPelhm=? “HEr PRINT USING “6k. fuse ype. eo a DO ele 
BE=? 
INPUT “ENTER QUTPUT VERSION (O=LONG.1=SHORT,2=NONE)", Tov 
BEEZ? 
INPUT “SELECT (O=LIO, 1f=VAP.2=(LIQ+VAPY/2)" [lay 
Dae vorgte 0 eo,e) TUBES 
ELECTRIC HEATED “ODE 
Di=Qrameter at thermocouple positions 
BAe Oa oe ee ae , 0129540, .072446, .0129540. O1guass 
DATA .0100965..01157..01157 
READ Jace) 
Di =) taCree) 
DZ=Diameter of test Section to the pase of ans 
DATA 015875, 0158752 et aee 5S, .015824..0156/5.9 015624. ee 
DATA .O127. ie. ro 
READ ace 
O1=inside diameter aunennancea ends 
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1925 Err2=tTv-Dla2z 

330" oleae = 1s 

1935 IF [nm=0 THEN PRINT USING “4%. 5¢4D0D 00.2%) 203 aia 
1940 IF [Rm=i AND I[de=0 THEN PRINT USING “5X. /CHDDSDD 22 0a 


UMNO. 1 1 Mbet wrot Le 
1945 IF Thm=! AND [dp=4 
Ale Pldee ty el sume er on 


THEN PRINT USING “2k SCH 0 aope 


PSS SAA ee 

FSS5e “GUiioe= 

1969 ELSE 

(965 IF ABS(At!ld-Dtld)>.! TREN 
1270" - deere oo 6 coe roe 

127 35> Joe or ecw 

1980: (BEe2 3000. -2 

IS8S° ELSE 

Peo) sBeee 300722 

1395 BEEP 30D0..2 

2000 END tr 

200S) (erel=atra-— iia 

eCid) Bidl=4ria 

eulS (e2r2Z='9- Jide 

220. -Dkazaty 

eUz25S IF ihm=) “HEN PRIST YSING “4x. SCNDD. OD lek) Oe ae ee 
2030 IF [hm=1 THEN PRINT USING “4X%.S5¢40D.00.2%) 22a D071 eee 
Penne. i milebs 1 ote, 

20>. (Meee 10 

204) GOTO ‘6385 

evas Ele 

2NSQ BEEP? 

COS5) Errl=Atia- Jig! 

“960 Olal=s4tia 

EU, sere = be ae 


68 


2.17. SUMP 
dnc tan 


ane “A ~ - 
1 2t.0c. ies 


aX) ,3(M3D.0D.2X) " ; Geis 


vs PSO 


T 
r r 
Lea eee ae 


= 
Lie 

2 
Oe iV 


2470 Jiazetv 

rag) ee tee fi Bs 1G 94% .S(400. 5D. 2K 2Dtie. tal. ide. iv. isu 
eve0 DF ihm=1 THEN PRINT USING °4%.3¢(ND90.90.2X%)":Otia.Fidt.Tldé.Tv.isum 
motile. lout 

2085 WAIT 2 

2990 GdoTUD [685 


meu eRRGR TRAP FOR Ido QUT OF BOUNDS 


IF Ido>2 THEN 
BEEP 
GOTO 1250 
Sip iF 
t 
fie eeeenin fF in-)) COOP 
IF Ikol=i THEN 2165 
Bee? 
iMUt CENTER BULK QIL 4” .Boo 


Tkol=1 

IF Ihm=) THEN OUTPUT 709:"AR SFO ALII VRS" 
IF Ihm=i THEN QUTPUT 709:"AR AFD ALS VRS” 
TIF Ihm=)) THEN Ntc=!2 

FOR I=1 TO Ntc 

GUTeuUtT Vos AS SA” 

Sum =) 

Be> «FOR Ji=! 10 20 

eeu ENTER 7O9:E 

2295 Sum=Sum+E 

2210 IF [='(9-Nsup) OR I[=¢t0-Nsub) THEN EtClirld=e 
met > NEXT Ji 

mee) Kdl=0 

Bees iF t=(9-Nsup) OR [=(10-Nsub) THEN 

2220 Eave=Sum/20 

B55 Sum=1). 

e241) FOR Jk=) TO 19 

2245 IF ABSCEt&( lk) -tave)<5,0E-5 THEM 

2250 = Sum=Sumrl ct (Uk) 

me 2 eee 

2250 Kdl="“dl+i 

Beo5> END IF 

2c7) NEXT Jk 

meee 8 =(2-—Ncun) JR [=(10-Nsup) THEN PRINT USING "4K%.°"Kdl = °" DD: Kdl 
e250 [IF Kdl>10 THE! 

Bees 3EEP 

Be30 BE=? 

ee eee Ge oA. 1oo mucn Scatcerina 1n data - repeat data set” 
Bevo 68GOT 1345 

meu> =ND IF 

eo END iF 

eS | Emt (TI) =Sum/(290-“d1l) 

Mee) NEXT I 

eoes )«€6UTF 6Ucfhm=! THEN 

233 OUTPUT 709:"AR AFZ9 ALCO VRS” 

Becs OUTPUT 709:"AS SA” 

2340 = Sum=0 

2345 SOR Kk=i TO 20 

Meo) ENTER 709:e2 

ecco Sum=Sumte 

eeou NEXT Kt 

2265) Emt (7) =8BS(Sum)/29 


LO COMONNM AION EP GW — 
DNGDNODMNODMNONeCWGMNO WI 
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= \Dilees, 


IF Inm=) “RE’, 

JUTPUT 7O9:"°3R APi2 AL‘ 3 VRS" 
rFoOR 22" -'8 2 

NUTPUT 709:"4S sa” 


Sum=) 

FOR Jizz! FO 2 

ENTER 709:E6 

Sum=Sumtreé 

Nie ee 

TF f=! THEN Vr=Sum/2 

IF IT=2 THEN Ir=Sum/2 

NEXT 1 

EHD 2F 

Pease 

IF [hm=0° THEN ENTER SF: le2:Son.foldS.Emr (<=) Yeoie 
ieee THEN “ENTER BFitee: 800, toidS.Emt (=? .-me 
END [ 


CONVERT emr’S TO TEMP.VOLT.CURRENT 
Twa=0 

FOR T=i TO Nte 

IF Iatc>) THEN 

LP -letatcl OR l=tbdter HEY 
TO) =233.99 

GOTO 2545 

END iF 

END TF 

TF Itt<4 AND Ihm=9 7HE! 

IF I>4 AND [<3 THEN 

CC =o oe oo 

GOTO 2545 

END IF 

END IF 

IX Pee N vs7tent<1)2 

NEXT I 

IF [t¢<4 THEN 

FOR f=" 7-4 

bP A =Wate | UR Vaslkatece Teen 
Twa=lwa 

Eoe 

Nia=iiat tC 1) 

END iP 

NEXT © 

Mia= las Co-rdte) 


TF fhm=! THEN 2650 
FOR f=t 10 3 
te! OR L=tdte2 THEN 


iva= Wat lt 1) 
evn) os 
NEX? 
Twa=Twas 3-Late? 
EMD TF 

Ld=' (9-Nsup) 
Tld2=TC10-Nsup? 
lidaa uid ld2) <5 
Ty=i (ii -tHeup? 


70 


Meee oe tk S SND Tame TRE 

Mo,5 23 la 2=- 3S, 23 

2530 Sere hee ct*))/2 

635 END Le 

25630 fF Thm=) THEN 2710 

e625 fsuno=T(12-Nsup? 

2700) Tainitet='Ct3-Nsup) 

2705 Toutz=t(14- Ee LD ? 

ead) iP ihm=) THEN 

2715 Ame=4BS( Tr) 

2720 Volt=ABS(Vr) «25 

2725) 6Q=Volt*Ame 

meow END IF 

2735 6oIF Itt=) AND Ihm=0 THEN 

2740) Keu=FNKcul Twa) 

Ges ELSE 

meee Keu=Kcua( ltt) 

eye END IF 

e750! 

2765? FOURIER CONDUCTIGN EQUATION WITH CONTACT RESISTANCE NEGLECTED 

2770 fF Inm=0 THEN tw=iwa-9-L0G(D2/D01)/(2-Pl+KcuKlL) 

Pi7> tk flav=0 THEN jtsat=!lda 

e7eo iF Llav=1 THEN Tsat=(Tldatfv)+.5 

Ome IF Llav=2 THEN Tsataitv 

27230 «iF Inm=1 THEN 

Ses Vava="inlet 

2300 Grag=37. 9853+.104388-1ava 

2380S Tdrop=ABS (Emr (7)) <1 .5+6/010"*Grad) 

Came) «=«6llavac=(inlet-Tdrop*-.5 

2315 IF ABS(Tava-Tavge)>.91 THEM 

See lavg=(tavarlavac)~.5 

See GUIDO 2800 

Zoeo END IF 

2350 ' 

2840' COMPUTE WATER PROPERTIES 

2845 Kw=FNKu( Tavq) 

2350 Muwa="NMuw' Tavg) 

205 Cew="NCow( Tava) 

2360 Pru=-NPryt Tavq) 

2865 Rhow="NRhow( Tavgq) 

fee lur=lavg 

oes ' 

<=380! Comoute MDOT 

=p Mdor=3.965/7E-3+Fms=+(3.919S5E-3-Fms-(3.3820066-56-Fms*' 1 .23688E-7-Fms+4.31397 
-19))) 

ea Seem Gon) Ucoo-iinlet*<1.36544E-2-‘' inlets. 252E-5))/1 .0037 

Kate! 

G20) wW=i*doteCou~tdrop 

See Patd=!droe/LOG((Tinlet-Tsat)/(Tiniet-Tdrop-!sat))? 

2910 VUo=Q/(PI-+Do-«L-t_mia) 

Bes Nw=D0*L0G6Do/D1)/(2. -Keu) 

Cee) iw=tsattFreLmea 

aa Yw=4dor/(RhoweP Ti), 2/4) 

eee! if Kett=0 THE? 

mea! Kat=! 

e941' Fdrop=idrop-.004"Vy 2 

eee GOD 3770 

Sesvee => IF 

2955 Rew=RhowseVy*401/Miwa 

Beem ete is/DisRew .2-P re (1/5.)-¢ Maas PE NMuw( Ti)) 214 

Soo mic='avq-}/(Pl-Dg-L-4.? 


a 


_—~e we 


J 


--+ -~4 —d ~—e =) ws —-O =e = bh 


OG) CIO ODIO CUI GO GOO CO COICO COCO CUIGOIC 
a ae 
NOomMonoulawMowowieo vo wo Ul 


SSO US On OO ~j Jone) Cll ee Co CoP) = 


Pesos ot wie 
Mgt Clete ie eee 
oo by 

END IF 

forbs urs loac) eS 

Ho=|/Ci/Uo-00/ (0144) =Ria? 
Thetap=}/ (Ho-P T<D)o-L) 

tw=lsattihetap 

Thetap=!w-tsat 

IF Thetab<0 THEN 

Beer ; . 
INPUT “TWALL<TSAT (O=CONTINUE, 1{=END)", lev 
IF Tev=3 THEN GOTO *32 

[F Tlev=!l THEN 3535 

END IF 

END IF 


COMPUTE VARIOUS FeOres I ies 
DY eee ae lsat) <=. 
Rhoze Rho film) 
Mu=FNMuCtrilm) 
K==NK(TFilm) 
Ce=FNCp(Tfilim) 
Bera="NBeta(Tfilm) 
Hf q=FNHfalTsatd 
N1=“4u/Rho 
Alpna=*/(Rho*#Cp) 
Pr="i/Alpna 

Peat ==NPsattisat? 


. | ' + 


mm 


m 


COMPUTE “IATURAL-CONVECTIVE dEAT-TRANSFER COEFF ECIENT 
FOR UNENHANCED ENDS) 

Hbar=!90 

Fe=¢(Hbare*P/(Kcu*9)) >. Seu 

Tanh=FNTanh(Fe) 

Theta=ihetap*!ann/Fe 

Xx2(9O.3!l+Setat*thetap*Do 3«Tanns(Fe«Ni-*Alpha)) .16 
Yo=Ci eC 559/P ro {S/o ows 7) 
dbarc2</Do-(.57.387**%x/Yy) 2 

[FP ABS CCHbar-Hbarc) /Hbare) > U0 Re IHE} 
Hbar=(4bartHbarc =. 5 

GOTH 2735 

Eno Le 


“J 


0) 
Oy 
> 


COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END<S) 
Ql=(Hbar<P<=KcureA) .S*Thetap*ianh 
Se-82 ei. 


ear Pe eck 


COMPUTE ACTUAL HEAT FLUK AND BOILING COEFFICIENT 
Qap=tc/As 

Htupe="Jaep/Tnetap 
Csr=(Cp-Therap/Hta)/(Odp/ (Muedfa) -¢ 014/79. 31eRho?) .5) 


RECORD TIME JF OATA TAKING 
IF m=) THEN 

QUTPUT 709:"™D" 

ENTER 709: 7ola$ 

CMD IF 


Te 


ee UT” arta 7G GEER 
ee PR EWTE~ 2S 7s 
Skat TF lows THE’ 
Sec ORIN, 
ee eon ix. Nata set adumoer = “"“.SDD.2X.°°8ulk Gi, 4 = ““.D0D.D.54, 5 
Bee .D0OP.1013¢ 
3295 iF ifhm=) THEN 
Bout) PRINT USENG why 1G. TOS ! 2 zg 4 . 5 ? 
SO TN TOUS eNG “10%. temp 3°", o¢?X. 400.90)": TO1), T02).103) 764). 705). 708), Te 
foe CO) 
Seegeees iit USING “10X."" Twa Tliad (vad. Seiya r © -sat lsum  : 
ees NG 10%-2(MDD.00,1X).1X%. MDD DD. 1xX.20Ci14.MDD.0D).2X.MDD.D°: twa. Tid. 
idee ty.esat.isump 
Soe PRINT WSING “10X.°° Thetab Htude Bada *™ 
Seem re iN! USING “10X.MDD.3D.1X.MZ.3DE.1X.NZ.3DE”;: Thetab.Htupe. lap 
SS Eis 
Eeogee PRINT USING “10X."" Fms Vu foac (inl drop Thetap q Uo 
ono. 
Peo ORIN) USING “*0X,4¢20.90D,! ote Dil AoC ODE, 1A) SeMS. Jia esate 
niet.idrop.ihetap.0¢d 
Sea>, END IF 
Sg50) END IF 
3395 IF [ovet THEM 
3360 IF J=1 THEN 
3365 PRINT 
3370 IF Thm=9 THEN 
Seo fl) USTHG “i0X%."" RUN ‘lo D1i24 Tsat 4tupe dp Thetao’’’°” 
3380 aS 
Seemeea iit USING “i0x,°" FMS ee Seat. HTUBE NDP THETAB’’* 
2330 END iF 
foo END IF 
3290 IF Thm=0 THE? 
MPMeree nemo ji iG) 1 274,00,4X.00.2h. 700.00. 3¢1X.4Z2.-SD0E)":J.800.isat.Htupe.Udp. Ine 
cap 
eae ELSE 
ems iicg 6 | 72k OD. AX. 0D.2X.400.90.3014.42Z.SDE) “iF as.30p.Jisat.dtupe.Jdp.f 
hertap 
Saou) =ND iF 
3425 END iF 
3430 fF im=) THEM 
Soee BE=>? 
Seem INPUT “OK TO STGRE THIS DATA SET (1=¥.0=N) 2° Ok 
Sees END [fF 
24S) IF Oke=i OR Im=! TRE) J=J+1 
3455 IF Ok=1 AND Im= 0 THEN 
S460 (fF thm=') THEN JUTPUT 2Filei:300.foiad.Emr (+) Vr, 
3465 IF fhm=1 THEN 3UTPUT Wrtriel:Bopo,ToidS.cmt (+) eas 
few END IF 
fom LF lure! THEI) GUTPUT sUF slesYu.Uo 
3480 IF [re=1 THEN DUTPUT aRer tle:ias. Ray 
Seog fr (im=! JR Ok=!) AND Iptotc=! THEN QUTIPUT 2Plot:Oqo. thetab 
3490 IF Im=0 THES 
eaee) 3EEP 
2300 INPUT “SILL THERE BE ANDTHER RUN (1=7,0=N) 2?" Go_on 
505 Nrun=J 
2510 teoolon=)) THEM 2535 
epee Go _on<>0 THEM Repeat 
Ee) Be 
Peewee oF J<‘lrunt! THE Reoeat 
2) END ES 


io 


-— - . —- 


Pe oolme CHE 
Sa2- 
aol SING “'OX. °° OTE: “",22.°" aata runs were stored am > tee 
tles 
aSSici ae (edi 
DUTP Sete Pano 
ASSIGN aFilei TO Ot_Fii 
ENTE Sei lei@Oatessudte! .lotece se. 
BUTPUT Spleens tes Ldtet cee 
FOR [=i TO Nrun-% 
[If Thm=93 THEN 
ENTER SF rlel:800.TolaS.Ear(«) Ve vir 
OUTPUT a elo soo lolat enn et 
PESe 
ENTER 3F rlet:Boo.TolaS.cmr(«) ,Fms 
QUTPUT @File2:Boo.TolaS.emf(*).Fms 
END IF 
NEXT I 
ASSIGN *File: TO * 
PURGE GUMMY" 
END IF 
BeSe 
PRT 
IF IplLot=1 THEN PRINT USING “72k.°°NOTE: “*,.2ZZ.°" X-Y pairs were storequaam 
dataurwle: “. 7 

ASSIGN ar: lez 13) sa 

ASSIGN aP lor TH = 
IF ITut=! TREN ASSIGN esUfirle TO < 
IF Ire=1 THEN ASSIGN sRerile TO «+ 
Cale Scats 
BEEP 

INPUT “LIKE TO PLOT DATA (1#Y,U=N2?~. dk 

IF Qk=1 THEN CALL Plot 
SUBEND 


5 


oy 


CURVE FITS OF PROPERTY FUNCTICNS 
DEF FNKeulT) 
DFHC COPPER 2590 TR 300 <x 


Wet trio S earn. 
K=434- .t12-ik 

RETURH « 

PES 


DEF FNMu¢T) 
(70.79 a <-CORUE St) Ve VS Cousi ey 


Tk =T+273,15 ‘CTH ¥ | 

Muze (80-4, 46567°¢19011.47/7Tk) 14! .OE-2 
Re Oe sail 

= MED 


DEF FNCo¢T) 

30 O40 Ue eres 

Lee ee oS aie K 

fp=,40188+! .6S007E-3-Tk+! .S1394E=551k 2-5. 58e5cE Ine ee 
Co=lo771 000 

RETURN Cp 


FNEND 
DEF FNRho<T) 
We 1*27e 25 ne Age 


Ra Cee tks fae oo A ry 8 

Ro=76. pea teal t~(°C1/2) +16. 413015-trt7, 476838"<>.5+!,) 1 a92ge4 
R9=909/ , 060425 

RE-URM Ro 


74 


o> 21) 
38275 
3341 
3345 
25) 
e255 
3360! 
3865 
3370 
oa/> 
3380 
3885 
go 20 
Bo95 
e200 
S205 
Bo 10 
6215 
B92 

e225 
3930 
e925 
3940 
e245 
6951) 
6255 
e260) 
$265 
eo 7) 
e275 
3989 
e265 
2390) 
e395 
4090 
4005 
4091 

4075 
429 
4025 
40350! 
4035 
eS) 
4045 
4y50 
4055 
46) 
4965 
4079 
4075 
4080 
408 

4999 
ay93s 
4100 
105 
4110 
aris 
i 


a 


=NEND 

DEF =NPec 7) 

Pee" NCoc Tyee NMuac Ty /ENK OT) 
RETURN 2r 

-NEND 

DEF FNK(T) 

ios0 < Wile | ING 
K=.071-.000261 +7 

RETURN X 

FNEND 

DEF SNTanh(X) 

PeeXP(X) 

Qe1/P ~ 

Tann=(P-0)/¢(P +0) 

RETURN Tanh 

FNEND 

DEF FNTIwsev (V) 

Pimorce, CC/) ical 

T=C(0) 
BOR [=1 TG 7 
Tat+C Ci) a) 7 
NEXT I 

toe tcal =i eihen 

T=1-§. 7422934E-2+7T+(9.0277043E-3-14(-9.3259917E-5)) 
BSE 

f=1+8.626897E-2+1T#(3./6199E-3-1+5.063892S9E-5) 

END ele 

RETURN T 

FNEND 

DEF FNBeta<T) 

Rop=FNRho( T+. 1) 

Rom=>NRho( T-.!) 
Betra=-2/(RoptRom)*(Rop-Rom)/.2 | 
RETURN Beta 

FNEND 

DEF FNHFg(T) 

Hf g=i .3741344E+5-7T-+(3.309436ic-2+Ts! 
RETURN 4fg 

FNEND 

DEF FNPsat( Te? ; 
feiwesteaego CURYE -il GF Peat 
aioe Gc Hoe 

Pa=5 249525771 re€.15552082+T Fr <0) W438 40963E-3+TF<9.5150671E-5)) 
Posra-i4./ 

IF Pq>0 THEN ' 
Psat=°4 

RECE 


Peat-—-g725.92/14./ 

END [IF 

RE“URN Psat 

FNEND 

Dec Insmootmn’ <,.s0p.isaz) 


vatoolas) 


‘Se oto. = lay Ag 


Wee o So). 005),005S).0¢9o) 

Pee OS 25 coo... 2 2048... 555e2.. 79909. (00252 
MA Pelee oe: oi oo ec eo os bees... 64197 
Denote 7725, .291 82. 94006, .31916.1.0845 
Se todee els... feso5..59609! ..605867..61889 
READ Ace) BO) Cle) DC) 

IF Bop<5 THEN L=Bop 

TF B3op=5 TRE!!! I[=4 


eo 


ee) 
ee 
ay ay 
414s, 
4150 
4155 
4150 
4165 


0495 
6500 
6505 


55402 


Dons 
6320 
6525 
6530 
6202 
6540 
$5545 
6550 
555) 
5555 
6550 
6565 
63/0 
bof2 
6360 
556 | 
pa385 
pe! 
Oo 
5600 
b605 
6610 
5615 


ce eat 
0O0<¢ 


S0c5 
6550 
boo) 


_—~- 


L= 3067} <2.) I=s 

ip isat=’ “den 
Hs=£XP(967)+89¢1)«L9G60%)) 
Bese 
dseekXP( CCl) +01) =1oG¢4)) 
END fF 

RETURN Hs 

FNEND 


SUB AL Ssonwc..G.) 

CEM Anlywenr Da DO. wu. eu 

OLM emt 12) 

ee PLOT SUBROUTINE GE es ie. ce earle . 
ieee 

INPIT "ENTER DATA FILE MAME”. Files 

BREED 

PRINTER IS { 


PRINPSUSING “ax, Siset “oom ontis 
PRINT USING "ak, He Q Vary Cs ana fecintaagee 
PRN SSG 4x." | “Fax (Oe Vary ere 
PRENTPUSING “4x7. “02 Vary Ce Pix 
PRN eUS LG ak os Cea Cl Pixel 1 


INEOIS ENTER OPT LUN fie 1x 
PRinNvter IS 70 
TF Tefitx=0 THEN 6585 
TF Icfix>0 TREN SEEP 
IF Icfix=1 THEN INPUT “ENTER CF" .CsF 
LP fesse THEM let “ENDea eer 
IF Icfix=3 THEN INPUT "ENTER CF, Ci’. CsF.Cu 
Paes boo 
INPUT wants io Vary Coot. 7. "=" J] = ceoer 
IF Iccoef=' THEN INPUT “ENTER COEFE’,R 
PRINTER IS 701 
Per x=) OR ee hee eT Cras .004 
betes p= 7 THEN CfazCs f 
Crawies 
Xn=.2 
Ef= 3 
Jy20 
RE So 
IF Tecoer=! THEN Rr=RX 
PRINTER IS 
PRIM Do.01.<cu 
ASSIGM @File 10 Files 
EN TEs wees be cNrin: Jaret cate | Bedtc.. —= 
rRw=Dor_UG(DosD1)9/¢ 2s" cud 
=1) 


aT f=! TO Nein 


SNTER BF pie:3on.TotasS.carle).Ems 
CONVERT EMF 'S TQ "EMPERATURE 

(ig =" les 

TOS) SENT vevCEmr6.l)) 

| Se a ees 

iSsat= Ce 2? nes 


Tava=T(S) 


76 


Meo) Grac=5/ 20537. -143559-' ava 

p74Q0)} Tarop2=Emr(7)4i 2 r6/ C70. *Graa) 

oye) ‘avac='(5)-farop-.> 

6750 IF ABSCiava-tavac)>.9)1 TREM 

eee iavg~' lavgriavgc)+.> 

eeoeu. GOD 67S 

6/65 END iF 

5770! 

5775! Compute sroperties of water 

§780 Kw=rNKw(Tavq) 

6785S Muwa=*NMuuw( Tava) 

6790 Cow=FNCow(Tavq) 

6795S) Pru=FNPru(Tavq? 

6800 Rhow=FNRhow( Tavg) 

§305S! 

5810! Compute properties of Freon-114 
poeemeemtda~ [droe/LO0G<¢(1(5)-Tsat)/¢(1(5)-Tdrop-Tsat?)) 
6620 IF Jy;=0 THEN 

6825 Tw=Tsat+Fr«imntad 

6830 Thetab=Tw-Tsat 

682 Jysi 

6840 END IF 

6845 Tf=(TutTsat)*.5 

6350 Rho=FNRho( Tf) 

6855 Mu=FNMu TF) 

63860 K=FNK(TF) 

6865 Cp=FNCp<(TF) 

6870 Beta=rNBetracTr) 

6875 Hfq=FNHfa(Tsat) 

6380 Ni="u/Rho 

6885 Alpna=*/CRhKo*C—? 

63890 Pr=Ni/Alpna 

6635! 

5900" Analysis 

poo COMPUTE MDOT 

Semoee A=" 1-(Do 2-91 °2)/4 

6915 P=PT«Do 

Gee Mdot=3.3637E-3+fms-(3.51955E-3-Fas7(8.32006E-5-Fms-(1 .23688E-7-Fms-4. 31397 
--10))) 

S225) 9=4dot*Cou-larop 

B220' COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
eee. FOR UNENHANCED SDS) 

63940 Hbar=!90 

53945 Fe=(Hbar<P/(Keu*4)) .5ebLu 

6350 Tann=FNTanh(Fe? 

5955 Theta=Thetap<«Tann/Fe 
6360 Xx=(9.31+Betarttherap*Do 3*tfann/ (Fee: <Alona)) .166067 
meee ly=(i+(.5939/Pr)°¢(3/%10)) (38/27) 

5970 Hbarc=</Dor(.5+.3587*%x/Yy) 2 

6975 IF ABS((Hbar-Hbarc)/Hbar)>.001 THEN 
6380 dHbar=‘(HbartHbarc)<+.5 

B25 GOI0 6945 

eoou «END iF 

Bess! 

7900! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENDS 
7005 O1=(Hbar<?=KxcurA) .5*Thetap*lann 
79190 %e=9-2-91 

Cimesenis=” +) et 

7020! COMPUTE ACTUAL HEAT FLUX 

7025 Qap=Ne/As 


ia 


woe 
ies 
7) 40 
7945 
7050 
0o5 
7160 
(G65 
7070 
7075 
7080 
ToS 
7090 
L0s5 
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OMmonononowm 
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as S) > LOUD 


Sod ee ae 


= NS ei NS 


_— 


CR Se- a) SR fee oe Sree 
Thetapslst /Cos4fa-(Qap/ (Mue)eds a) = 
Ho=30p/ hetap 

Omeqa=4o/CFf 
Uo=N/(P1+Do-L-tatd) 
Vuw=4dot/(Rhow*P 1<+D1 2/74) 
Reuw=howsrVy*) 1 /Miuwa 
Twir2=lur)*Ru/ (Pl «Dowd 
Gama=K/DieRew .3*P rw €1/2.)<(Muwa/FNMuw(Tw1))>.14 
PRINTER) 1S 5) 

Yy=(1,/Uo-Ru) * mega 

Xw=Jmeaa*Do/ (Gama*D1) 

Sx=SxrXu 

Sy=Sy+rYu 

Sxy=Sxy+Yu <u 

Sx2=Sx2+Xu+Xw 

Sy2=Sy2+Yur ly 

NEXT I 

ASSIGN aFile TO * 
M=(Sx#Sy-Nrun*Sxy)/(Sx«Sx-Nrun*Sx2)?) 

C=(Sy-Sx=) /Nrun 

IF Icfix4#0 OR Ictix=3 OR Icfix=4 THEM 

Creai79 

Cfie=17C 

END IF 

Ir Lefix=! THEN 

Cie=1/" 

Cre=c* 

env te 

if Gerry =2 “HEN 

Crem 

Che=i1/7 

END IF 

Ee tetre=-2 (HEN 7260 

[TF ABSC(C;-Cic)/Gic)>.001 UR ABSCCC! -Cted Cred > 00 ae 
Caet Ores, 5 

Crate +t fee Ss 


~~ 


So? ee 
O° 4/09. @terho)) (50 ae 


PRINTER IS ¢ 
PRINT USING “OX."" Caf = °".MZ.3DE.0%."" Ci = "42 SDE eee 


PRINTER 25 3701 
GOTi si35 


END TF 

Peay 

PRINTER 1S.701 

PR WSchGreecc kee (Cues 

PRIN USING. 3X. HHACSUMED MZ. SDE .o%a"7.20E4 0 eee 
POINT USING Hay "eral CULATED AL CDE. GAs eee ae 
PRINT 

oum2 =Sy2—2M=Sxy- 2-0 SS y oN Stoke Sos ee 

PRINT USING “10X.°°Sam 6+ Squares = "G7, 20e = sme 

PRINT USING “0%. “Caer? tetent. =". 40 4000s a. 

SUBEND 

oe a 

Ameo (le one 

Mie 2. AE -S-10)° 4 

RETURN “uy 

= MEND 

DE -FNCow(7) 

Cow=4,21120853-<(2. 260258 -3- 1904 4206 1b See oe eee 
fe WRN Cowes 1000 

-“IEND 
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moo) 
7340 
7345 
iS) 
P35 
7360 
7505 
SAD, 
"S75 
7230 
7365 
73390 


PUCOS 
moe 710 
D375 
1230 
125 
0390 
ross 
10400 
10405 
10410 
10415 
10420 
MnyaeS 
10430 
10435 
10440 
19445 
10450 
no455 
10460 
1465 
10470 
e475 


DEF ="IRhoaw( 1) 

Rm ee se a ee ea toes | CE =- 3-71. 234147E-5)> 
RETURN Ro 

- NEND 

DE= SNPerwT) 

Pru="*ICowt) =F NMuw( T)/EN Ku CT) 

RETURN Pry 

FNEND 

DEF FNKw< 7) 

Keres. our ee o..05 

hue oe o/h at 283 I5-<*( i e00/-X<¢ .52577- .07384*X))) 
RETURN Kw 

EMenO 


SUB Seiect 


COM / Ldpy Tdo 

SEEP 

PRINTER IS + 

Solis USiNiGwee eso rect Operon: 

we iGo. | hakeina Gata or re-processing erevious data” 
Pili Uotlomone | Plotting data on Log-Log ~”” 
Bo momeiicemowe ee: Plotting data on Linear = 
Puli wistiiemegk seo take cross-olot coefft file’’” 
Pell] Msp on. “4 Re-cireculate water’ 
FRibieuoums OA. 5 Purge” 

Reh Use SX2°°5 f=Deooecorrection’”” 

PRINT USING "6X%.""7 Print Uo File” 

INPUT dp 

IF fcp=0 THEN CALL Main 

IF Iqgo=!1 THEN CALL Plot 

Me teo=2 THEM CALL Plin 

TF tap=3 THEN CALL Coer 

TF Tap=4 THEN CALL ifain 

iF elapes TmeEN CALL Purg 

ie Soap=5 THE CALL Tdaen 

teeeiao=, (HEH CALL Uoprt 

SUBEND 
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APPENDIX B 


FLOWMETER CALIBRATION 


A Fischer Porter Flowmeter was uSed in the experimental 
apparatus to indicate the water flowrate. Prior ee 
conducting the boiling tube data runs, calibration of the 
flowmeter was performed over a temperature range from 19 °C 
to 38 °C, The goals of the calibration were, 


1) to develop a mathematical expression relating flowm- 
eter percent reading to mass flowrate and 


2) to determine a viscosity correction factor due to 
varying temperature range and flowrates. 


A weigh tank, platform-type scale, and stopwatch were used 
to record water collection data. 

Data were processed using the program "FMCAL", listed 
later in this appendix. This program took the flowmeter 
percent reading, water weight collected, and elapsed water 
collection time as the inputs. The outputs were an experi- 
mental mass flowrate andamass flowrate difference. This 
difference was given by a mass flowrate computed from a 
correlation minus the experimental mass flow rate for the 
Same flowmeter percent reading. The correlation was based on 
a fourth-order least-squares fit to the calibration date 
After all the data runs were completed, the data points were 
combined into one file and a single fourth-order polynomial 
was generated to describe the relationship between the 
flowmeter percent reading and the mass flowrate. Of the 
original data points, 98.61 percent were within a 9235 
percent range. Reprocessing the data points led to /78 
percent within a + 1.02 percent range. The equation (B.1), 
resulting from the calibration data analysis, is a fourth- 
order polynomial from a curve fit of the reprocessed data 


without a viscosity correction factor: 
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oo eco es 04S) | (3.1) 
4 
-8  S020CMMES)- + 1.2369 
=) 51 OG eeenneS) 


(FMS =flowmeter 
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Experimental discharge rates are less than the rates 
calculated from theoretical flow equations primarily because 
Pemeete E@rtect of internal molecular friction of the fluid 


Paemmete Viscogity. The theoretical flow equation is: 





ecV .(@%-0 ) a 
8 eer — > - | ea 

ie eS 4 
Boaeme O=volumetric flowrate, A =cross-sectional area ofctne 
Marrowest part of annulus, Y-=volume of float, f@-= density 
ee Loat, Ppl -™eos-seclilonal area of the largest pazxt of 
moat , and “9, density of San. Mine V\Gcetcicient of 
Peel = 22, Cy mS a Sumerewomeof Viscosity. A constant value 


for C im equation (3.2) would be desirable as this would 
Smowea Tlegdation in the variations due to the viscosity. The 
Seeeee- Of the float shape on C, with varying Reynolds number 


? 
ay 


XS 
MemcmOwl im Figure 3.1 (a). The Reynolds number for the 


flowmeter uses an effective diameter of the tube inner diam- 
S@ememinus the float outer diameter. This plumb bob style 
Meee presents a large surface area to the fluid stream and 
MueMmesCOUSs €fleccs increase as the flow rate increases. 
wees Jl: pp. 9804-9808] 

Mime) data colleéted from the calibration runs showed 
Ppeete or no change in the Discharge Coefficient over the 20 
Bem) pezcent flowmeter setting range and th 15 .cKk venance 
me 2luid temperature. Daeiecguarce Sedge plums bob float 


MeereS da Constant Coefficient af Discharge at Reynolds 
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number > 10000. This is due to the effective shape of the 
float on the streamline pattern as shown in Figures B.2 (a) 
amd (b). The experimental calibration range of 1300 < Re < 
27000 displayed a flattening out of the Coefficient of 
Discharge. A different effective shape is seen by _ the 
streamline pattern because of the minimal clearances between 
the tube bead guides and the float. Also, the float used in 
this particular flowmeter had tapered edges versus square 
edges. This different effective shape is shown Figures B.2 
(c) and (d). Figure B.1(b) shows the effect of the float 
shape on the Coefficient of Discharge for the tapered edge 
EVoat . 

Additionally, the accuracy of the flowmeter contributes 
to the flattening effect on the Coefficient of Discharge 
shown by the experimental data. The accuracy effect is felt 
in two major components, reproducibility and scale factor. 
Reproducibility accuracy is minimally determined by the 
change in the fluid flowrate, corresponding float movement, 
and observation of float position. The calibration data had 
good reproducibility over the various flowrate and tempera- 
ture ranges. The scale factor accuracy is inherent to the 
scale markings on the flow tube. Fischer Porter conducted 
numerous tests of flowmeter calibration readings. They used 
a maximum scale fraction error of 0.9 mm on a 250 mm tube. 
Based on their results, Fischer Porter reported [Ref. 3l: 
pp. 9814-9816] an accuracy of plus or minus one percent at 
high flow rates was a reasonable expectation. Also, at low 
flow rates (< 20 percent) a resonable expectation would be 
an accuracy of eight to ten perecent. 

The experimental data followed a trend of negligible 
viscosity effect. The flow tube design, float shape, and 
scale factor accuracy were the main contributors to this 
Meedelon of a yiscosity correction factor. As such, equation 
(Bel) was used without a viscosity correction factor to 


calculate the mass flowrate for the boiling tube runs. 
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eee bay 
Boba e2ie we esON PLOT 


AseEeterencsa im @hapters Ill and IV, this appendix 
provides the data reduction procedures used to computes 
imside and outSide heat-transfer coefficients. The tubes 
Maecaerotescddta colleceion were the smooth, modified Hizh 
Peete. , and Modified TuxsSo-3 tubes. The modified Wilson plot 
anaylsis could not be directiy run on the enhanced tubes 
Beemee crelzjable boiling corztelations for these surfaces do 
mot exist. Also, only two unknowns can be allowed between 
Seerelations for the inside and outside heat-transfer coef- 
ficients. The steps used in this procedeure are outlined 
Below and aze contained in DRP4S, - which is cited in Appendix 
A? 


mas c Equations 


ietemoverall. hedt-transter coetfticzent in terms of the 


Meesall €nesmal resistance is given by: 
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mF oe Rt Cou) 
O° Paes i. "O 
where 
es een a = (C. 2) 


MOC ol. ed Wilson Plot 
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1. Assume Coe and Te¢ for the boiling side and ecoemeuee 








wes L] afoy-9,,) - oo (C.35) 
: (2: | e091 4 nS. 
ae a 
2 
or = (C. 4) 
O (c_.)t ; 


Note: with r equal to 0.333, equataen bay? represents the 
Rohsenow correlation {|Ref. 23:7 @emme7G 7). 


2. Assume Cj; and compute: 





a.i4 
a. a oo) peeeeereo Lu 
ae D Re ie CC. 5 
ad Ww 
ne = Co] (Com) 


J a COUR StisuGe on and hj from steps 1 and 2 above @iimea 
equation cc. 1) and rearrange to yueld. 








= sy 
Sa (co)? (Cc. 
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(C88 
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Construct the least-squares line for Y versus X neem 
On mo ©: 


Y=mxX +c (Cc. aoa 
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Peas eo blows: 


Da Somouce a mew Sseceerevalues for C5 and C_ 
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6. Repeat steps 1 through 5 until convergence for Cy and 
Ceof¢ between two successive lterations is less than 0.1 
percent. 
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APPENDIX D 
DATA RUNS 


The following table outlines the data runs used in this 


investigation: 
Where: 
i, WH - Smooth tube 


in 


oO CO N OO 


LO. 


HF- High Flux (i.e. , porous-coated Korodense) tube 


HFM - Modified High Flux (i.e. , porous coating machined 
off) tube 


TB - Turbo-B tube 


TBM - Modified Turbo-B (i.e., external enhancement 
machined off) tube 


T - inlet water temperature 


V - water velocity 


const - constant 
dec - decreasing 
ine - increasing 
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Description 
Poor en beet 
ee eee 

dehe Off Effect 
warm water 

1eht Off Eitfect 
7 eam ) 

const V dec 

a0 const V dec 
ar const V dec 
Ay const V dec 
i aout V cee 
ah const V de 
Repeatability WHO 
A eons & : de 


de 
eae OFF REESE 
warm water ) 
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ik const V_ dec 
C2 ste Off Effect 
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cold wat 
ight Off E 
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epeatability 
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const V 
const V 
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Repeatabil 
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Const V 
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dec V 
Ime V 
Leder-Tate 
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leder-Tate 
leder-Tate 
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Coeff 
Coeff 
poene 
eff 
Reece 


Feet 
fect 


of HF12 


HF 28 ) 


Run 
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Const Vv 


const 
const V_ dec 
inc V const 
const V dec 
dec V const 
const V_ dec 
ine V comer 
const V dée 
const V 


const 
const V_ dec 
inc V comst 
const V dec 
dec V const 
inc V const 


const V 


Hn H- 
ae 


cold wate 
LeneeOr t E 


a ada ad aa ae 


leder-Tate 
Sieder-Tate 
Sieder-Tate 
sLleder-Tate 
Sleder-Tate 
Sieder-Tate 
Sieder-Tate 


warm water) 


dec 


he eee Effect 
nee fete Eo 


fect 


Coeff 
OSes 
Coeff 
Coeff 
Coeff 
Coeff 
Coeff 


de 
Repep tabi ey Se TBOsy 


dec 
eee of TBIZa 


Run 
Run 
Run 
Run 
Run 
Run 
Rive 


eae Uncertainity “for mass flow rate, Reynolds number, 
heat flux, er, ~Wwal resistance, OVverahewnede-coansier 
@eemricient, inside heat-transcter coefficient, and outside 
heat-transfer coefficient were analyzed for selected runs of 
mee smOOch, High Flux and Turbo-3 tubes. The analysis was 
based on the Kline-McClintock [Ref. 32: p. 3] method of 
Meicestainity analysis. For example, the following equations 


were used for the uncertainity of the heat flux: 


Q 





and 


Memeo corcancs Wien Aline and McClintock, the umcertainitzies 
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. 
genuine asrccainicy in the canduction losses for the unenhanced 
Peea Of tae Boiling tube were considered negligible in 
Pemeaemsom ©<S sae bolling surface wuncertainity and as such, 
Reeem@Giesnezareaed. “alable 5 lists the uncertainities for the 


previously mentioned values. 
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UNCERG ADIN? 7: eo wipveeeoees 
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fue WHOS WHO 5 HF 13 foes T304 TB04 
ae 31140 §©6©26690 49540 38780 63780  -§3670 
w/m 
Deis 1.39 eee 1. 39 1. 387 1.07 
(m/s) 
= L. 56 oes L. 56 2.5m 1.43 2.5m 
sRe a ; 
=e 1.79 tees Lars 2uen tevone 2a 
5 
= 1.59 25/0) eso ee 17a 22 
4796) | Sais, 10 eee 7 DS lL. 66 
sR 
= 2.76 2.76 fa wee? 9 eel omy 2.92 2.92 
ao) 5.21 . 4.39 3.48 ‘eee al 3 G2 
ne 
= 2s 28S evoe 2S L. 48 .. 2a 
rhs 
L 
a 6945 5.08 9 4.90 Jaman 3.35 3.96 
A 
= 2 
Note: All scuns were verlormed in pure Setsetzecame 
Na. = Smootawtupe 
He = Korodense tube wth Higa nu Coa ean 
[3 3.31200 -2 eo oe 


The large uncertaintty Eor the wall fesvsecamee Bo 


due to theme 


HF13  (verous- cosaced KXorodense tube) is 
percent umcertainity of the thermal concuctivic7 Sco ee 
cient. For this investigation, the coefZicient Used eae 


average value for the values given in the 


[= 
open literature 
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ARID E: 
biot OF NOMENCLATURE 


NOMENCLATURE 

A - Surface Area 

CoH Specific heat at constant pressure 
ee ot Concentration 

Manco lae oteder-lTate-type coefficient 
wep pg nonsenow coefficient 

d; - Inside diameter of tube 

d, - Diameter of mouth of reentrant cavity 
D - Diameter 

e - Depth of internal ridge or fin 

Be ene rey transfer per unit volume 

foes Touctilon factor 

zg - Acceleration due to gravity 

<a Gzeavitational constant 

h - Convective heat-transfer coefficient 
i Seeecir lc enthalpy of vapor-liguid mixture 
k ¢ Conductive heat-transfer coefficient 
ky - Conductive heat-transfer coefficient of metal 
L - Tube length 

le = EMNtrance effect length 

meee- tiass flow rate | 

p - Distance between ridge or fin peaks 
fe - Prandtl number 

me - Heat flux 

eee total heat transfer 

mee - lube radius 

Re - Reynolds number 

[wees omermal wall resistance 

wee otanton number 
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U, - Overall heat-transfer coerivervense 
Vz - Specific volume of vapor 
V - Velocity 
= (Chiaiwee (ieee T . ae eee 
GO - Surface tenSion 
Oe qe oUt Lace= bens 1 ommo fcael 
Gg. > surtace tCensilonvot tren meer ams 
Uo =" Viscosity. 
Ua) Vs COS tiymG Loaner 
Ht - Viscosity of refrigeramt 
9 - Density 
Pace Densiey ot Jo14. 
a> Density of retrigerame: 
SUBSEREPTS 
fee= st lid 
1 - inside 
in - inlet 
m - mixture 
Oe Oues Tae 
@ l=. Oleh eau id 
out- outlet 
rl- refrigerant-liquid 
S - Smooth tube 
Sdb-WSeLurarvion 
w - water 
wo - outside wall 
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